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We study the spin-boson model with a sub-Ohmic bath using a variational method. The transition 
from coherent dynamics to incoherent tunneling is found to be abrupt as a function of the coupling 
strength a and to exist for any power < s < 1, where the bath coupling is described by J(u>) ~ 
auj 3 . We find non-monotonic temperature dependence of the two- level gap K and a re-entrance 
regime close to the transition due to non-adiabatic low-frequency bath modes. Differences between 
thermodynamic and dynamic conditions for the transition as well as the limitations of the simplified 
bath description are discussed. 



I. INTRODUCTION 

The spin-boson modeli is a paradigm model for the 
study of dissipation and decoherence in quantum me- 
chanics, and as such it is has been applied in a wide 
range of systems. Such applications include the search 
for macroscopic quantum coherence^, electron transfer 
in chemical^ and biological^ physics, and most recently, 
the problem of dephasing and relaxation in solid state 
qubits&. 

The particular case of the sub-Ohmic spin-boson model 
has had an interesting development in the last few years. 
Compared to the Ohmic bath, the sub-Ohmic bath is 
characterised by an increased density of states for the 
low frequency bath modes. This makes analysis of the 
dynamics difficult as the low frequency modes generally 
lead to non-Markovian dynamics and strong memory ef- 
fects, even for relatively weak coupling. One of the physi- 
cal situations corresponding to the sub-Ohmic bath is the 
1// noise in Josephson qubits^. Although there are cer- 
tain limitations and assumptions in describing 1// noise 
by an equilibrium sub-Ohmic batb-4, the study of the 
sub-Ohmic spin-boson model may be useful in this and 
other contexts. 

In the earliest treatments of the sub-Ohmic modeli, it 
was argued that the sub-Ohmic bath always destroys the 
coherence of superposition states and localises the system 
in one state for any non-zero coupling. This conclusion 
was based on the non-interacting-blip approximation^ 
which fails in the weak-coupling limit to the bath. More 
recently, several works^^S addressed the problem of the 
sub-Ohmic spin-boson model and found that coherent 
phases can exist for sufficiently weak coupling. 

In the light of these developments, we contribute to 
this discussion of the coherent dynamics of the sub-ohmic 
model by demonstrating the existence of the coherent 
regime for arbitrary s < 1 using a simple and intuitive 
variational method. The variational method we use was 
originally developed by Silbey and HarrisiS for the prob- 
lem of Ohmic damping. We will show that such a treat- 



ment allows us to define precise criteria for the thermo- 
dynamic existence of a coherent phase, and also provides 
a means to quantitatively map out the parameter space 
of the sub-Ohmic coherent regime. Within the coher- 
ent regime, we also give new results for the renormalisa- 
tion of the parameters that describe the coherent dynam- 
ics. A re-entrance regime close to the coherent-incoherent 
transition is found. In addition, strong coupling to non- 
adiabatic modes is considered, showing the limitations of 
Silbey-Harris variational ansatz. 

In section II Al we briefly give an outline of the spin- 
boson model and in section^ we describe the variational 
method and explain the simple physical picture behind 
the variational ansatz. Throughout this paper we are pri- 
marily concerned with finding the conditions under which 
coherent oscillations of the TLS arc possible, and in sec- 
tions lIIll and llVl we give some new quantitative results for 
the critical couplings and renormalised parameters of the 
dissipative tunneling at zero and finite temperature. In 
sections IVI and IVI1 we highlight some of the limitations of 
this method and in section lYlIl we discuss our results and 
compare the findings with the results obtained by other 
authors. We end with a brief conclusion and summary. 



A. The spin-boson model 



The spin-boson model consists of a single two-level sys- 
tem (TLS) coupled linearly to an infinite bath of har- 
monic oscillators. The TLS can be thought of as span- 
ning the two lowest levels of a double well potential, or 
in other contexts, the TLS may appear in the situation 
where the transition matrix element between two given 
energy levels is much larger than the transition matrix el- 
ements to all other energy levels of the system. The two 
levels are coupled by a tunneling matrix element K, and 
taking these levels to be eigenstates of a z , the spin-boson 
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Hamiltonian is given by, 



H = Kd-r 



i 



uiffidi + 1/2) + d z ^ 9i{di* ■ 



-ai) 



(1) 

e is the bias energy between the minima of the wells 
and for the rest of this paper, we set e = 0. dj,ar are 
the bosonic annihilation/creation operations for the bath 
modes. The last term in EqQ, the linear coupling of d z 
to the coordinate displacement of the bath oscillators, is 
assumed. 

In the absence of coupling, the eigenstates are coherent 
superpositions of the left and right states, and the parti- 
cle can oscillate in time between the wells at a frequency 
2K. The couplings {<#} between the spin (or TLS) and 
the oscillators generally lead to damping of this motion, 
or may even suppress tunneling entirely. For an in-depth 
discussion of the rich dynamics of spin-boson models, the 
reader is referred to the original review of Leggett et ali 
or the more recent collection of papers on this subjecliii. 
The observed dynamical behaviour is determined by the 
spectral function of the bath, J (to), 



J(uj) =n^2gfS(u 



(2) 



For frequencies below a high energy cut-off, u) c , the spec- 
tral function can be modeled by the power law form, 



(3) 



where a is a dimensionlcss parameter that measures the 
effective strength of the system-bath coupling and u> s is 
an energy scaled included to keep a dimensionless. In 
this paper we assume that the cut-off frequency oj c is 
much greater than all other scales in the problem, but 
the scale w s can be large or small compared to other 
scales e.g. K, KbT. 

The case of a bath with s = 1 is known as the Ohmic 
bath, and the dynamics and thermodynamics of this 
model have been studied extensively in the literatueAi 
Baths described by s > 1 are termed super-Ohmic and 
we will not discuss them further. Here we shall focus 
on the sub-Ohmic spin-boson model where the bath is 
characterised by < s < 1. 



If the tunneling is switched back on we propose that 
the approximate eigenstates of the system correspond to 
a dressed particle tunneling between the two levels car- 
rying a dynamical cloud of such oscillator displacements. 
Our variational ansatz is therefore in the spirit of an adia- 
batic approximation^ but with a special treatment of the 
low frequency, non adiabatic modes as we discus further 
in section IVTT1 

To describe this physical picture, we reproduce below 
the technique originally used by Silbey and Harris^ for 
the ohmic bath. We begin by performing a unitary trans- 
formation of the spin-boson Hamiltonian Eq.Q 



H = UHU' 1 
The unitary operator U is given by, 



U = exp 



-d z Y^^t 1 fi ( d i _ ^ 



(4) 



(5) 



The arbitrary coupling parameters /; introduced in 
Eq.lJSJ are proportional to the effective displacement or 
dressing of each bath mode due to the coupling to the 
TLS. If fi = gi then the transformation diagonalises the 
last three terms of Eq. , but as we demonstrate in sec- 
tion i|Vll|) . setting /; = gi is often a sub-optimal choice 
for {/*}• 

As we mentioned in the introduction, we are primarily 
interested in establishing whether or not coherent oscil- 
lations can exist in the sub-Ohmic model and to answer 
this question, we introduce the quantity K which is given 
by, 



K = K ( cxp{ 



i 
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bath 



(6) 



where the brackets denote the thermal expectation 
value taken over the bath modes. 

We interpret K as the effective coherent tunneling ma- 
trix element of the dressed particle. The exponential fac- 
tor in Eq.© suppresses the bare tunneling, and this fac- 
tor arises due to the partial overlap of the oscillators that 
dress the TLS as it tunnels between the wells. Adding 
and subtracting K to Eq.QJ, we re- write the Hamilto- 
nian as 



II. THE VARIATIONAL METHOD 

In this section, we motivate the variational approach 
to this problemiS and outline the method. In the absence 
of the tunneling term, the spin-boson Hamiltonian Eq.JQ) 
reduces to the well-known independent boson modeliS, 
and the solutions of Eq.QJ correspond to the particle lo- 
calised in one of the wells. The oscillator part of the 
Hamiltonian is then just a collection of displaced oscilla- 
tors and can be diagonaliscd by a simple translation of 
the oscillators oc gi(d z ). 



H = H Q + V (7) 
Where the separation is into a main coherent part Hq, 

H = kd x + J2u l {d l U l + 1/2) (8) 
I 

+ Y,(fi -2/, Si), 
i 

and a series of perturbation terms, V, which contain the 
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remaining weak coupling between the TLS and the bath, 



V 
Va 



V+ 



V+a 



V— U- 



V <?z, 



(9) 
(10) 



Xcxp{-2^ fi^idi - d; f )} - £(11) 



The introduction of K and the separation of the Hamil- 
tonian into a coherent part and perturbations is remi- 
niscent of a mean-field type theory, and our treatment 
is essentially in this spirit. Considering the main part 
of the Hamiltonian, Eq.JHJ, we see that if K is finite, 
the eigenstates of are coherent superpositions of the two 
levels and the TLS can undergo coherent oscillations be- 
tween it's levels. If K vanishes, the degenerate levels 
become uncoupled and no coherent oscillations are pos- 
sible. Therefore, at the mean field level, we can use the 
existence of a finite effective tunneling matrix element as 
the signature for the existence of a thermodynamic co- 
herent phase in the TLS. The point where K vanishes as 
a function of the system parameters marks the transition 
to the incoherent phase. 

However, this thermodynamic criterion for distinguish- 
ing between coherent and incoherent dynamics is only 
approximate, as we have not yet included the dynamical 
effects of the perturbation terms on the TLS dynamics. 
The effect of the perturbations and the general limita- 
tions of using K as the criteria for the transition from 
coherence to incoherence will be discussed in section (|Vj) ■ 

To calculate K, we must first determine the {/;}. Fol- 
lowing Silbey and HarrisiS, we compute the Bogoliubov- 
Feynmani^ upper bound on the free energy of the system, 
Agi^- Bogoliubov's theorem^ states that the true free 
energy, A, of the Hamiltonian J7J) is related to Ab by, 

A < An 



A B = -/T 1 lnTrcxp(-/?F ) + (V) Ho + 0({V 2 ) Ha )- 

(12) 

Due to our choice of i?o, we have constructed the per- 
turbations so that {Vi)H = 0. In section lVll wc explicitly 
calculate the second order terms and find that they give 
a small contribution to Ab if oj c — ► oo. Therefore, drop- 
ping all higher-order terms, Ab is given by, 



A B = -K B T\n[2cosh(KP)}+J2^ 1 (f?-2f l g l ), (13) 



where we left out a term due to the free bath ground state 
energy which does not depend on {/(}. As can then be 
minimised by varying the {/;} to find, 



fi = gi [1 + 2KUJ- 1 coth(w| / 9/2) tanh/3X) . (14) 



Notice already at this stage the limiting behaviour of 
coefficients {/;}: 



fl 



gi if uji/3 > 1 and K(3 <C 1 
gij~ if ^iP > 1 and lui <C K 



(15) 



The coefficients of effective coupling {/;} vanish in the 
limit loi — > and finite K . 

We now substitute this form for {/;} back into Eq.® 
and use the spectral function, Eq.@, to turn the sum 
over modes in Eq.JU) into an integral. We then obtain 
our key equation, the self-consistent equation for K 



F[k] 



K = Kexp(-2F[K\) 



J{uj) coth(^)dw 
2 J RTtanh(/3 K) coth(^)) 



(16) 



(17) 



III. RESULTS AT T = 

The appearance of K on both sides of Eq. (|16fl means 
that we must solve self-consistently for K. At T = 0, 
we can perform the integral in Ea. 1)16(1 exactly. The re- 
sults below were calculated by extending the upper limit 
in Eq. (|17fl to infinity, an approximation that be easily 
dropped but which is a possible approximation for s < 1 . 
Calculating the integral and substituting into Eg. 11611 . 
we obtain, 



K exp 



((2A > ) 1 " s sin( 7 rs) 



K. 



(18) 



The self-consistent values of K can then be obtained 
by numerically solving Ea. p8|l for general values of a. 
Note that K = is always a solution of Eq. (|18|) . 

Using Ea. ((T%|l . it is possible to determine the criti- 
cal coupling strength a c for fixed K. a c is the cou- 
pling strength above which the only possible solution of 
Eq.® is K = 0. 

To see the existence of this critical coupling, we define 
the LHS of Eq.lPSJ as <f>(a,K). This function has the 
typical form shown in Fig^ and crucially, has only one 
minimum for any sub-Ohmic bath. Finite solutions for K 
exist when the <j)(a, K) intersects with the line K = K , 
and the point of intersection is controlled by the coupling 
strength a as shown in Fig.Q. The critical coupling 
strength can then be clearly identified as the coupling 
strength where the minimum of <p(a, K) just touches the 
line K = K as shown in Fig.QJ. When a > a c there is 
no longer any intersection with the line K = K and the 
only self-consistent value of the rcnormaliscd tunneling 
matrix clement is K = 0. 
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FIG. 1: Solutions of the self-consistent equation for s = 1/2 
and K/u! 3 — 100. Solutions exist where 4>{K) intersects the 
line K = K. As the coupling a is increased, the curve shifts 
until there is no intersection. The coupling where this occurs 
defines the critical coupling a c . 



The position and value of the minimum in (f>(K) can 
be determined by elementary calculus, and this gives the 
results, 



2K n 



a c A(s)Y /(1 - S) 



(19) 
(20) 



<f>(K, 



mint — J-^min 



where A(s) = eir s(l — s)/ sin(7rs). From these equations 
we then find that the critical coupling a c is given by, 



1 / 2K 
A{sj \~ 



(21) 



Note that unlike the case of ohmic damping, the sub- 
Ohmic critical coupling depends on the ratio of 2K/u> s , 
and that the Silbey-Harris approach predicts a finite a c 
as s — > 0. Note also that the above condition fEal21|l can 
be rewritten so it is a condition on the coefficient aw\~ a 
appearing in the spectral function (Eq|3J). 
We also find that when a < a c , K satisfies, 



K exp{- 



1 



1 



} < K < K 



if 



a < a r 



(22) 



This inequality shows that at T = 0, K undergoes a 
discontinuous jump from a finite K to K = as a —* a c 
Only at s = 1 does this method predict that K ^ 
continuously. 

In some other treatments of this problem^, the energy 
scale u> s is set equal to the high frequency cut-off uj c . If 
this is done in this method, we find that wc cannot send 
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FIG. 2: Behaviour of K(a) as a function of a at T — 0. For 



this computation K = 10 



and lo c — * do. 



u> c — v oo in Ea. l|17|l as this leads to K = for all s and 
a. Keeping tu c = tu s finite, we get the result, 



A(s) 



2K 



(23) 



The only modification to the previous result, Ea. (|21|l . 
is the exponential factor and the replacement cu s — > cu c . 
The exponential factor is the correction for a finite cut- 
off and is important only as s — > 1. Note that as s — > 1, 
Ecr. (12*31) correctly predicts the critical coupling for the 
ohmic case, a c = li We will not be too interested in 
the ohmic case as this has already been thoroughly dealt 
with in the literature. Therefore for the rest of the this 
paper we work with s < 1, u> s ^ lo c and oj c — > oo in the 
integral Ea. ("T7Jl . 

For a ^ 1 we can make a perturbativc expansion in a 
and determine K{a) to first order in a, 



K(a) = K 



1 



a sir r u> s 
sin(7rs) \2K 



(24) 



For general couplings to the bath, equation ((TTj|l needs 
to be solved numerically for K, and the typical behaviour 
of K(a) across the whole range of coupling strengths is 
shown in Fig.®. 



IV. FINITE TEMPERATURES 

A. High temperatures 

Calculating the integral in Eq. i|17|) and solving Eq. 1161) 
for the general case of finite temperatures can only be 
done numerically. However some analytical results can 



5 



be extracted in certain limits. For the case of high tem- 
peratures and very weak coupling where K j3 <C 1, we 
again find a coherent regime which crosses over to inco- 
herent relaxation at a critical temperature T*, 



K 



«/(*) 



2K 

us* 



a < a c 



(25) 



where f(s) is a slowly varying function of s which is al- 
ways « 0(1). In this regime, we find that the transition 
from finite K to K = occurs discontinuously at T*. 

For stronger coupling the relation given by Eq. l|25|> is 
violated, and a numerical study we have performed shows 
that K vanishes at a significantly lower temperature than 
T* as a — > a c . For weak coupling, the numerical calcula- 
tions of K(T) give values of T* in good agreement with 
Eq.JSJ. 



B. Low temperatures 

For temperatures close to zero where Kf3 3> 1, we 
can solve the self-consistent Ea.l|16|l for weak coupling 
by making a perturbation expansion in powers of a. The 
result to first order is, 



K{T,a) = K(0,a) + 2ag(s) 



(KbT) 2 ( 



(26) 

where K(0) is given by Ea. (|24|l and g(s) is another func- 
tion of s which is of order unity. Ea. H26(l shows the 
surprising result that K becomes larger as the temper- 
ature is increased from zero. This result was also de- 
rived by Weiss^ for the ohmic bath and was qualita- 
tively described by Kehrcin and Mielke^ for the sub- 
ohmic bath. However we believe the quantitative result 
given in Ea. (|26|l have not been explicitly presented be- 
fore for the sub-ohmic bath. We shall discuss this effect 
in more detail in section IVIII 




FIG. 3: K as a function of temperature for a sub-ohmic bath 
with s — 1/2 and a range of different couplings. For this 
numerical computation K = 10~ 3 o; s and uj c — > oo Note that 
the two lowest curves have finite values of K only between 
T* and a lower re-entrance temperature. For these curves, 
a > a c w 0.021. 
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FIG. 4: Schematic plot of the boundary between coherent and 
incoherent regimes as found by the variational method. The 
curve is proportional to 1/a for small a. The plot also shows 
the small re-entrant region for coupling strengths greater than 



C. Intermediate temperatures 



For intermediate values of a and T, we can deter- 
mine K numerically and the typical behaviour is shown 
in Fig.©. In all cases we find that K increases to a 
maximum and then drops discontinuously to K = at 
T*. In section IV 111 we estimate that the peak in K(T) 
should occur approximately at a temperature KBT max ~ 
K{T max ), which for sufficiently weak coupling can be ap- 
proximated as KBT max ~ K . Comparing to the numeri- 
cal results we find that this is a good order of magnitude 
estimate, but the peak typically occurs at a lower tem- 
perature ~ T max /2 as illustrated in Fig.©. 

The numerical results also reveal an interesting feature 
if we look at the temperature dependence of K for sys- 
tems with a > a c . We find that for couplings slightly 
above a c , the TLS is incoherent at T = 0, but then 



develops a finite K between some re-entrance tempera- 
ture and T*. In this re-entrance regime, K shows the 
same non-monotonic temperature dependence described 
above, and some examples of the behaviour of K in this 
"super-critical regime" are shown in Fig.{jy|. 

As the coupling between the TLS and bath is increased, 
the re-entrance temperatures and T* merge to one finite 
temperature and beyond this coupling, K = for all 
temperatures. This region is generally very small and 
together with the results for a < a c , wc obtain the 
schematic coherent and incoherent regions of the sub- 
ohmic model as shown in Fig.Q). This re-entrance phe- 
nomena is a consequence of the same mechanism that 
causes the enhancement of K at low temperatures, and 
wc discuss this effect in section IVIII 
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V. TLS DYNAMICS AND LIMITATIONS OF 
THE VARIATIONAL METHOD 

In section [H] wc defined the criteria for coherent dy- 
namics as the existence of a finite renormalised tunneling 
matrix element K. However this criterion does not take 
into account the effect of the perturbation terms given 
in equation These perturbation terms introduce 

dissipative dynamical effects which can alter the oscil- 
latory behaviour of the TLS in the coherent tunneling 
state. These effects can be calculated by a variety of 
methodsiiSiLiSi and for the ohmic case are well under- 
stood. 

For the sub-ohmic problem, we are interested in the 
weak-coupling behaviour where approximations like the 
non-interacting blip modeU^ are no longer valid. The 
weak-coupling should however permit us to analyse the 
effects of perturbations using the perturbative reduced 
density matrix method^i^ To second order in the per- 
turbations, the reduced density matrix of the TLS, 
p s (f),is given by, 



Ps{t) 



dt Tih 



V(t), V(t),p s (t)p b (0) 



(27) 



where the operators are written in the interaction rep- 
resentation V(t) — exp(iHot)V exp(— iHgt), and V and 
H$ are defined in Ea. (|llf) . pb(0) is the thermal density 
matrix for the unperturbed bath modes. Once p s (t) is 
known, all the observables of the TLS can be found us- 
ing (o) = Tr [p s (t)0], where the trace is only over the 
states of the TLS. 

As can be seen in equation H27|) , the time development 
of the reduced density matrix depends on the whole his- 
tory of it's motion, and such memory effects can lead to 
strong modification of the tunneling dynamics^. As a 
simple example of the dynamical effects that perturba- 
tion can cause, we consider very weak coupling and ignore 
the memory structure of the bath. This simplification is 
known as the Born-Markov approximtioni^, and apply- 
ing it to the spin-boson model, we find that a z obeys the 
simple equation of motion^*i£ 



dt 2 



2T 



d(d z {t)) 

dt 



+ 4K 2 (<f z (t)) =0. (28) 



Therefore, in the Born-Markov approximation, the co- 
herent oscillations of the TLS are exponentially damped 
with a decay rate given by, 



r, 



J(2K) coth(K /3) 



(29) 



so that as t — > oo the TLS settles into a decoherent mix- 
ture of localised states. The coherence of the initial state 
is gradually destroyed by interactions with the environ- 
ment on a time scale 1/r. This is a generic phenom- 
ena for open quantum systems^ and for long enough 
times, the initial coherence of the superposition state is 



destroyed. Therefore when wc talk about the coherent 
phase in the Silbey -Harris variational method, we mean 
that the initial ground state is coherent; the subsequent 
tunneling is then subject to decoherent and dissipative 
processes which eventually destroy the coherence. 

The purpose of these remarks on dynamics is to point 
out that in the thermodynamic coherent phase, these de- 
coherent and dissipative processes can potentially drive 
the coherent tunneling of the TLS to become incoher- 
ent. Therefore it is possible that there is a transition 
to incoherent motion due to dynamical effects that may 
occur before or after the thermodynamic transition we 
have found at a c or T*. For example, the Silbcy-Harris 
variational method predicts a c = 1 for the ohmic bath, 
whilst it is well known that for ohmic baths at T = 0, 
there is localisation for a > 1, incoherent tunneling for 
0.5 < a < 1, and coherent oscillations are only observed 
for a < 0.5i 

In order to find such a dynamical cross-over in the 
variational method, we need to account for the pertur- 
bation terms and solve for the dynamics of the TLS. For 
Ohmic damping, Silbey and Harris calculated that the 
crossover to incoherent tunneling occurs when the equa- 
tion of motion l|28|l becomes overdamped, which occurs 
at a coupling strength cti nc = 2/tt at T = 0. If we apply 
the same procedure to the sub-Ohmic bath then we find 
that, 



2K 



s-l 



(30) 



after substituting Eq[2sl|we obtaion the relation between 
dynamic «m C and thermodynamic a c critical couplings 



Ctinc S a _ 2es(l - s) 



sin(7rs) 



T = 0. 



(31) 



This result shows that a.i nc /a c ~ O(l) for < s < 1 
but different in general. 

However, this result only applies if the use of the Born- 
Markov approximation is valid and, except possibly at 
extremely weak coupling^ this is not a good approxima- 
tion for sub-Ohmic baths. The Born-Markov approxima- 
tion fails for sub-Ohmic baths due to the presence of low 
frequency modes which cause large correlation times and 
strong memory effects. It is generally recognised that 
non-Markovian effects lead to stronger decoherence than 
that described by the simple Markov rate, and at present 
there is much discussion of non-Markovian dynamics in 
the context of qubit decoherence ratesiSiiLi2iiS. 

As we mentioned in section^ our method is based on 
a simple and intuitive variational ground-state, and we 
have ignored the dynamical effects of the perturbations 
in our discussion of the transition between coherent and 
incoherent phases of this ground state. The existence of 
a finite K as the signature for the coherent phase can 
be thought of as a thermodynamic criterion for the co- 
herent phase, and in light of the discussion above, the 
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critical couplings we have deduced from thermodynam- 
ical considerations can be different from those deduced 
from dynamics. 



VI. CORRECTIONS AND FAILURES OF THE 
VARIATIONAL GROUND-STATE 

The determination of K relies on the minimisation of 
the free energy bound, Ab, given in Eq.|0. In this 
section we estimate the higher-order corrections to this 
free energy. We have already shown that the first order 
term in powers of the perturbation vanishes and so the 
first corrections are given by second-order terms. The 
second order term in the Bogoliubov-Fcynman bound on 
the free energy^ is given by, 



.4 



(2) 




> WH °Ve- WH °VdW 



(32) 



The calculation of the second order contribution to the 
free energy is outlined in appendix^ For weak coupling 
at T = 0, we find that the contribution to the free energy 

is small, -g— ~ 0[(^-) s ] for uj c — > oo. Therefore we 
expect that our calculations based on the minimisation 
of Ab to be accurate in the weak coupling regime. For 
stronger coupling the corrections have to be calculated 
numerically, and again we find that corrections to Ab 
are small when ui c is much larger than all other energy 
scales. 

Another potential weakness of the method is that the 
variational ansatz may not be a particularly good guess 
at the true ground state in the first place. We can in 
fact demonstrate some cases where the variational solu- 
tion is sub-optimal. For simplicity, we shall show this by 
considering a spin-boson Hamiltonian with only one bath 
mode. 

We call bosonic modes adiabatic if the frequency of 
such modes is much larger than TLS frequency ojb 3> K, 
because these modes can follow the TLS adiabatically. 
The Silbcy-Harris approach is accurate in treating these 
adiabatic modes as well as being exact in the K = 
localized state. Now we turn to the opposite situation, 
the anti-adiabatic case, where K 3> £j&. 

We introduce a different variational wavefunction for 
the TLS in the basis of d z . It is given by, 



|*) = 



(33) 



where the number is a real variational parameter to 
be determined and <f> = corresponds to ||). Notice that 
in this ansatz we allow parity symmetry (up and down 
direction for the spin) to be broken unlike in the Silbcy- 
Harris approach. Wc fix the TLS in the variational state, 



and this gives us an effective Hamiltonian for the bath 
mode2£ given by, 



H, 



eff 



7— I? — 7a ) + +ua)a (34) 



This is an example of an independent boson Hamilto- 
nian and can be diagonaliscd exactly^. The resultant 
ground state energy is given by, 



E„ 



2K<t> 
l + <t> 2 



<£_ 

Id 



1 



(35) 



and wc minimise this energy with respect to <j> to find the 
optimal ground state wavefunction. The result is that the 
optimal value of </> is given by, 




(36) 



This result shows that the ground state spin is a linear 
combination of the localised and delocalised states c.f. 
the variational method where the spin ground state is 
a purely delocalised or localised state. The part of the 
wavefunction corresponding to the localised state gains 
a displacement energy, whilst the tunneling energy is re- 
duced as the tunneling part of the wavefunction has a 
reduced weight due to the normalisation of the wave- 
function. 

We notice that the important parameter here is 
g 2 /{ujK). When g 2 /(cuK) ;§> 1, what we can call the 
strong coupling case, the parity breaking (i.e. <j) + 3> 1) 
is large. Since we assumed K ^> u>, the strong coupling 
case implies that g 3> u). We will now show that when 
this strong coupling condition is met, the Silbey-Harris 
method gives a sub-optimal groundstate. 

For 2g 2 /loK 3> 1, which corresponds to either strong 
coupling or a very low frequency bath mode, the ground 
state energy of Ea. i|35|) is, 



E„ 



(37) 



The corresponding bound found using the Silbey-Harris 
method at T = is, 



-4, 



-K exp 



2K 2 



<r_ 

K 



(38) 



There is a large region of parameters that the ansatz 
of EqESI has lower ground state energy than the Silbey- 
Harris ansatz. In particular, if we set g 2 /K 2 <C 1, then 
if u> is sufficiently small so that w g, we find that, 



A j 



-K - %- » E g . s . 
K 



(39) 



Therefore, when these conditions arc satisfied, the Silbey- 
Harris variational method is sub-optimal. For constant 
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coupling, g, we always enter this breakdown regime as 
us — *■ 0. However, since the coupling constant g{us) can 
be frequency dependent, there can be (and are) many 
situations when weak coupling is valid asw — > 0, provided 
g(us) vanishes quickly enough to maintain g(uj) <C u>. 

These results show that the coherent state found by 
the Silbey-Harris method can be sub-optimal for baths 
with finite couplings between the TLS and low frequency 
modes. In appendix [B] we highlight this by comparing 
the variational ground state given by l|33[) and the Silbey- 
Harris state in the limit s — ► 0. 

VII. DISCUSSION 

In section [H] we stated that the physical picture be- 
hind Silbey-Harris approach is that the tunneling par- 
ticle drags along a cloud of displaced oscillators as it 
tunnels between the wells. For modes with frequencies 
much larger than the tunneling frequency we expect this 
adiabatic approximation to work well. The complications 
arise in this problem due to the presence of low frequency 
modes in the bath, especially in the sub-ohmic problem. 
These non-adiabatic modes cannot follow the tunneling 
motion and need to be treated separately from the adia- 
batic modes. 

If we try and treat all modes with the same adiabatic 
approximation and set /; = <?/, then it can be seen that 
the integral in Ea. H17fl diverges in the infra-red and al- 
ways leads to K = i.e. no coherent oscillations. This 
complete suppression of tunneling for the sub-ohmic bath 
was also obtained by Leggett et al* using the technique 
of adiabatic renormalisation. 

However, the variational method goes beyond the adi- 
abatic approximation and finds solutions with finite K. 
The appearance of a finite K can be traced back to 
the free energy bound we calculated in Eci. (|12[l and it 
is shown explicitly in Eq. i|A2(il) . There arc two compet- 
ing processes, the choice /; = g\ maximises the second 
term, the dressing/displacement energy. However, for 
sub-ohmic baths this always renormaliscs K to zero and 
thus incurs an energy penalty. Ea. (|A26|l is a non-linear 
function of a, K, T and which process dominates depends 
sensitively on these parameters. When a < a c (T) it is 
energetically favourable to have a finite K. 

For a < a c and T = 0, we see from Eq. l|14|) that the 
variational method has loosely separated the bath modes 
into two distinct sets. Modes with us > 2K respond adi- 
abatically to the tunneling motion i.e. have /; ~ gi- 
Non-adiabatic modes with us < 2K couple more weakly 
to the TLS, with coupling strength /; fts gij~ as usi — * 0. 

This vanishing of the coupling at low frequencies pre- 
vents the infra-red divergence in Ea. (|17f) by fixing an 
effective cut-off at 2K tanh(if/3). In this method, the 
free energy minimisation naturally determines the cut-off 
for the mode elimination, unlike in the adiabatic renor- 
malisation scheme^. We also note that while the non- 
adiabatic modes decouple from dressing the particle, they 



have not disappeared; they give the dominant contribu- 
tion to the perturbation term V$, Ea. HlUf) . and can cause 
significant dynamical effects. 

The variational method also predicts interesting be- 
haviour for K(T) at low temperatures. As we demon- 
strated in section llVl K(T) initially increases with tem- 
perature and this behaviour can be seen to arise from 
the non-adiabatic modes. Interestingly, we find that the 
renormalisation of K{T) due to the non-adiabatic modes 
actually decreases at finite temperatures. We would nor- 
mally expect that as the temperature is increased, the 
occupation of low frequency oscillators would increase, 
and this should lead to increased renormalisation through 
the hyperbolic cotangent factor in Eq. l|17|) . However, 
from Eq. fl^jl we see that the dressing due to modes with 
fcgT < us < 2K decreases with temperature, and this 
decoupling leads to an overall reduction in the renormal- 
isation of K due to these non-adiabatic modes. 

The dressing parameters for the adiabatic modes are 
effectively independent of temperature, and so when they 
are thermally excited they always renormalise K towards 
zero. At low temperatures when there are almost no adi- 
abatic modes excited, the reduction in the renormalisa- 
tion due to non-adiabatic modes leads to the increase of 
K with temperature. This low temperature reduction in 
the renormalisation of K also gives a natural explanation 
for the re-entrance of finite K at finite temperatures for 
systems with a > a c - At high enough temperatures, the 
renormalisation due to excited adiabatic modes always 
dominates and K then decreases until it goes discontin- 
uously to zero at T*. 

The exact point at which the adiabatic modes halt the 
increase in K depends sensitively on the relative weight of 
adiabatic and non-adiabatic modes and thus depends on 
the spectrum of the bath. However, we can still estimate 
where the maximum occurs. From the discussion above, 
the turning point occurs at the temperature at which 
the adiabatic modes begin to be excited. This occurs 
approximately at a temperature KsTmax ~ K(T max ). 

There have been several other recent treatments of the 
sub-ohmic problem^SiS and we find that this simple vari- 
ational method is consistent with several of the main re- 
sults. The flow equation analysis of Kehrein and Miclkc 8 
also showed that a coherent phase exists for the sub- 
ohmic model. On the basis of the well-known connection 
between spin-boson model and Ising model in statistical 
mechanics^, the coherent phase, corresponding to the 
high-temperature disordered phase of the Ising model, is 
expected to exist. Many results of Refi£ are in fact con- 
sistent with ours, including a qualitative prediction of the 
rise in K(T) at low temperatures and the discontinuous 
transition at zero temperature. 

It is important to remember that the transition in the 
spin-boson model as a function of coupling constant a can 
be related to the transition in an infinite one-dimensional 
Ising model with long-range interactions as a function 
of temperature, only when T — in the spin-boson 
model. Therefore comparison of the nature of the transi- 



9 



tion (lst-order or 2nd-order type) is limited to T = OK. 
In this paper we are mostly concerned with the transi- 
tion of spin-boson model at finite temperature, with sev- 
eral parameters describing the bath (a, u) s , cj c ). Yet the 
comparison with the results known for the Ising model 
with l/r 1+s interactions indicates that higher-order cor- 
rections to Silbcy-Harris ansatz should be necessary to 
describe the close proximity of the transition, since for 
s > the transition in the Ising model is of 2nd-orderi2i 
The numerical renormalisation group analysis by 
Bulla, Tong and Vojta- found that the system is localized 
at s = 0, and their pcrturbative RG results suggest that 
for s > 0, the transition is continuous as a function of a. 
As our method is based on a variational ansatz, we can- 
not make any strong statement about the exact nature 
of the transition. As we noticed in sectior lTTll there are 
several parameters which describe the bath, and the tran- 
sition may depend on the constraints between parameters 
imposed and assumptions used in the mappings to other 
models. 

Shnirman, Makhlin and Schorl have also demon- 
strated that coherent oscillations are possible in the sub- 
ohmic model, but their work focusses on calculating the 
dephasing and relaxation times of the dynamics rather 
than renormalisation effects. In contrast to Bulla, Tong 
and Vojta, their diagrammatic approach predicts that 
the TLS can be coherent at T = and s = 0. As we 
discussed, differences between thermodynamic and dy- 
namic properties are expected for the spin-boson model 
with a sub-Ohmic bath, and further understanding of 
these questions is desirable. 



VIII. CONCLUSIONS 



dependence of the effective dressing parameters I f/) (ll4|l . 
and we have highlighted the natural separation in this 
method of adiabatic modes (w > 2K) and non-adiabatic 
modes (u> < 2K) . Our numerical study of this theory also 
found a new phenomena, a re-entrant coherent phase that 
exists at finite temperatures for systems with a > a c if 
a is sufficiently close to the critical coupling. 

Importantly, we showed that dynamical and thermo- 
dynamic criteria for the transition are different and sen- 
sitive to non-adiabatic modes. We also discussed several 
limitations of the description of the spin-boson model by 
an equilibrium bath characterized by the spectral func- 
tion J(uj). 

We would like to thank P. B. Littlewood for useful 
discussions. 



APPENDIX A: CALCULATION OF 
SECOND-ORDER TERMS FOR THE 
FREE-ENERGY BOUND 



In section I VII we discussed the size of contributions 
to the free energy from higher order terms in equation 
(|12p. In this Appendix we outline the calculation of the 
lowest order correction terms to the free energy bound. 
The fist corrections are second-order in the perturbations 
and given are by, 



.4 



(2) 




,WH a 



Ve 



-WH 



VdW 



(Al) 



We have studied the sub-ohmic spin boson model using 
the intuitive variational method of Silbey and Harris^ 
This method has allowed us to reproduce a number of 
previously known results about the coherent sub-Ohmic 
model, but without having to make lengthy or unduly 
complicated calculations. With this in mind, we note 
that this method may be useful for a first look at different 
types of environment for which there is some question 
about the existence of a coherent phase. 

For the T = Sub-ohmic spin boson model, we have 
shown that coherent oscillations exist if a is below a crit- 
ical coupling, a c which we have explicitly calculated in 
Ea. (|2ip . When this condition is met, the renormalised 
tunneling matrix element K satisfies, Ke 7 ^ < K < K 
and undergoes a discontinuous transition to K = as 
a — > a c . 

We have also presented new numerical results which 
show the dependence of K(T, a) on temperature and cou- 
pling strength. We have shown that K{T) has a non- 
trivial dependence on temperature, initially rising to a 
maximum value and then decreasing to a discontinuous 
transition at a critical temperature T*. We were able 
to show that this behaviour arises from the temperature 



The perturbation terms are shown in equation and 
the Hamiltonian Hq is defined in equation (JHJ). The av- 
erage is explicitly given by, 



Tr exp(-(3H )A 
Trexp(-/3ff ) ' 



(A2) 



Each perturbation term is a product of a spin opera- 
tor and a bath operator. As the thermal density matrix 
corresponding to Ho is also separable into spin and bath 
parts, we can calculate each term in equation HAljl as the 
product, 



dW 



e WH o V s e~ WH ° V s 



e WH ° V b e~ WH ° V b 



(A3) 

here s refers to the spin part of Hq and b is the bath 
part. Before discussing these factors, it is useful to re- 
write the perturbations in terms of the spin components 
x,y,z instead of the raising and lower operators. This 
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*ives the perturbations as, 

Vo = (si - fi) ($i + tii ] )a z 



Vi = K 



cosh{2 ^2 fl ^i 1 (til - a V)} - 1 



(A4) 



Ox (A5) 



Vz = — iK s'mh 



2^2 fiui\tit-aj) 



o v (A6) 



1. Spin part 

The spin factor is of the general form, 



pjk = ( e WK<T i(T . e -WK, i(Tk } (A7) 



where i = x,y,z. The exponentiated spins can be 
written^, 



cxp(6<Ti) = cosh(#) + sinh(0)<7j 



(A8) 



and using this and the pauli spin algebra, one can derive 
the general relationship, 



= cosh(28)(jj + ieijk sinh(26 , )(Tfe 



— 2 sinh (O)Sij <jj . 
Substituting this into l|A7[) we get, 

P s jk = cosh(2K W){a j( r k/S 



(A9) 



+ ie ij ismk(2KW)(oiak)s 

- 2smh 2 (KW)S lj (a J a k ) s (A10) 

and finally, all the spin factors can be calculated using, 



{o x ) s = -tanh(^/3) 

(<Jy) S = (<T Z ) S = 0. 



(All) 
(A12) 



2. Bath factors 

If we define the operator ai(W) = 
exp(W Hi,)ai exp(—W Hi>), then the bath terms contain 
only averages of the form, 



(v l (w)v J ) l 



(A13) 



where the Vi are the bath parts of the perturbation terms 
defined in equation (|A6|) . To continue we need to calcu- 
late these expectation values. For the terms involving 
products of Vi,2 the following theorem is very useful. If 
the operators A and B are linear in the co-ordinates or 
momenta of an oscillator, then it can be shown^i 



(e A e B ) b = e h[( A2 )f+(B 2 ) b +2{AB) b ] 



(A14) 



for example, if we define A(W) = 2J2 l fiuJ l 1 (di(W) 
tii (w)), then 



(v 2 v 2 ) b = 4(r L ^ 



4 



A(WO)( e A(0) _ e -A(0)- 



2 

= -K 2 smh(-f(W)) 
where j(W) is given by, 



7(W) = (A(VF)A(O)} 

= -4^/,V[e-^ 



e (A(W)A(0)) _ e -(A(W0A(0)> 



(A15) 



(A16) 



-UJ t W 



(ni + 1)] 



and we have used K = K exp(i(A(0) 2 )). 

Combining these results with the spin factors, we cal- 
culate that the second-order contribution to the free en- 
ergy is, 



2A 



(2) 



- 

K 2 I d W [cosk{~i{W)) - 1] 
Jo 

■P f 

K l I dW sinh(70T)) cosh(2_ftT W) - sinh(2_ftT W) tanh(K/3) 



P 



dW{ V 2 (W)V ) sinh(2^ W) - cosh(2A' W) tanh(^/3) 



fP r . 

/ dW(V (W)V 2 ) smh(2KW)-cosh(2KW)tanh{Kp) 
Jo L 

rP , 
J2(gi - fi) 2 / dW (e WuJ, hi + (ni + l)e~ Wu}l ) (cosh(2KW) - tanh K f3 smh(2K W) 

i J o V 



(A17) 
(A18) 
(A19) 
(A20) 
(A21) 
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(2) 

Note that the free energy correction, A B ' , is stated for 

(2) 

the case of finite K. For K — 0, A B is identically zero. 



3. Second-order terms at T = 

At T = we can calculate all expectation values in 
(|A21I> explicitly and the free energy correction takes the 
form, 



APPENDIX B: ALTERNATIVE VARIATIONAL 
TREATMENT OF SPIN-BOSON PROBLEM 

In this section we give a treatment of the sub-Ohmic 
spin-boson model using the alternative variational solu- 
tion given in section IVII In the anti-adiabatic or non- 
adiabatic situation of K ^ 10, TLS can be thought as 
creating effective potential for bosonic mode. As before 
we write a variational state for the spin, 



.4 



(2) 



2 

2 



dW [cosh(j(W)) - 1} (A22) 
dW smh{ 1 (W))e- 2kw (A23) 



1 



1 



(Bl) 



o 



We then calculate (\&| H s b to get the effective Hamil- 
tonian for the bath modes. This is given by, 



- 6K 



; E 



{oJi + 2Kf 



(A24) 



H, 



2Ko 



1 



I 



,t, 



We will show that the typical size of the correction 
term is small compared to the main free energy in the 
limit of large u> c , which is the normal situation in this 
model. For the term (|A24J| we get a contribution of, 



gi{a + a f ) + u>i a' 
(B2) 

The first tunneling term is minimized for real </>, so that 
<f> is chosen to be real although in general complex. 

Again, this is a set of independent boson Hamiltonians 
and the energy of the variational ground state is given 
by 



9i 



{ui + 2Kf 



2K 



1 



1 + s 



s-2, 
(A25) 



E, 



g.s. 



2K(f> 
l + d> 2 



1 - 



1 



2 2 

I 



(B3) 



where we have introduced the spectral function and 
approximately calculated the integral. The other two 
terms. I|A22|) and (|A23l) . cannot be evaluated in a simple 
analytical form, but we note that as lu c —> oo these terms 
give finite contributions if s < 1. 

The main part of the free energy Ab is given by, 



The sum over the bath couplings can be explicitly calcu- 
lated by substituting the spectral function into the sum 
to get, 



E 



U)l 



aivt 



U3 



A B = -K + Y,{fi-1fi9l) 



1 — S f-UJe 



-K - 



(cj + AK)uj s duj 



o (w + 2Kf 



(A26) 



(A27) 



"27 



(B4) 
(B5) 



where again, we have used the spectral function to con- 
vert the sum into an integral. Under assumption K <C 
u) c , the leading term in lo c of Ab is 



A 1 



2s V u s 



(A28) 



(2) 



Comparing this to the second order correction Aw , we 
see that corrections due to the term given by l|A24(l are 

small, and are controlled by the small parameter {^-) s 
for s > 0. As we let oj c become large, the corrections 
from terms (IA23I and 1A22|I tend to a finite value, whilst 
Ab grows as uj s c . Therefore, the relative correction from 
(|A23|) and l|A22|> becomes small in this limit. However, 
these and higher order perturbations may still be relevant 
in the proximity of the coherent-incoherent transition as 
K -► 0. 



Looking at the ground state energy i|B3() we see that 
as s — > 0, the static displacement energy of the oscilla- 
tors (given by the second term of Eq.jBjjJ)) diverges and 
becomes the dominant term for any non-zero coupling. 
Minimising the free energy w.r.t <j> we always find that 
<t> = (or <j> = oo) and therefore the particle is always 
localised for any non-zero coupling at s — 0. This is due 
to the fact that these soft modes have no resistance to 
the static force due to the spin in the limit u>i — * 0. 

For s = 0, the Silbey-Harris variational method pre- 
dicts a coherent phase with finite K for sufficiently weak 
coupling. The free energy of this state is, 



-4, 



-K 



2 



-K — auj s In 



(uj + AK)duj 
(w + 2K) 2 

OJr. 



2K 



(B6) 
(B7) 



Comparing the energy of this coherent ground state to 
the energy of the localised ground state, we see that for 
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s = and uj c = oo the localised state is lower in energy 
for any non-zero coupling between the bath and the TLS. 
The coherent state is therefore never favourable when 
s = and the finite K found by the variational method 
is an artefact of the method. This artefact occurs due 
to the divergence of the static displacement energy of 
the oscillators (singular limit for lo c —> oo), which causes 
problems with the free energy minimisation we use to 
determine K etc. Notice though that K — is also a 



solution of the self-consistent Ea. (|17fl . and so the Silbey- 
Harris method can correctly describe the s = state if 
we ignore the sub-optimal solution with K > 0. 

To summarise, the divergence of the static displace- 
ment energy of the oscillators for s < implies local- 
ization in the ground state and dramatic differences be- 
tween thermodynamic and dynamic properties. Such dif- 
ferences due to non-adiabatic modes can also be seen for 
s > 0. 
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